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Kinetic theory of Dirac fermions is stndied within the matrix valued differential forms method. It is 
based on the symplectic form derived by employing the semiclassical wave packet build of the positive 
energy solutions of the Dirac equation. A satisfactory definition of the distribution matrix elements 
imposes to work in the basis where the helicity is diagonal which is also needed to attain the massless 
limit. We show that the kinematic Thomas precession correction can be studied straightforwardly 
within this approach. It contributes on an equal footing with the Berry gauge fields. In fact in 
equations of motion it eliminates the terms arising from the Berry gauge fields. 


1 Introduction 

Dirac equation which describes massive spin-1/2 particles possesses either positive or negative energy 
solutions described by 4-dimensional spinors. However, to furnish a well defined one particle inter¬ 
pretation, instead of employing both type of solutions, a wave packet build of only positive energy 
plane wave solutions should be preferred. A nonrelativistic semiclassical dynamics can be obtained 
employing this wave packet. Semiclassical limit may be useful to have a better understanding of some 
quantum mechanical phenomena. In the massless limit Dirac equation leads to two copies of Weyl 
particles which possess opposite chirality. Recently, the chiral semiclassical kinetic theory has been 
formulated to embrace the anomalies due to the external electromagnetic fields in 3 -|- 1 dimensions 
mm- This remarkable result was extended to any even dimensional space-time by making use of 
differential forms in [3] by introducing some classical variables corresponding to spin. Although in 
[3] non-Abelian anomalies have been incorporated into the particle currents, the solutions of phase 
space velocities in terms of phase space variables were missing. In [4j a complete description of the 
chiral semiclassical kinetic theory in any even dimension was established by introducing a symplectic 
two-form which is a matrix labeled with “spin indices”, without introducing any classical variable cor¬ 
responding to spin degrees of freedom. In this formalism, although the classical phase space variables 
are the ordinary ones, the velocities arising from the equations of motion are matrix valued. It has 
been shown in [5] that this matrix valued symplectic two form can be derived within the semiclassical 
wave packet formalism mm- In fact the “spin indices” label the linearly independent positive energy 
solutions. 

In the formulation of the Hamiltonian dynamics starting with a first order Lagrangian, the related 
Hamiltonian should be provided. In the development of the chiral kinetic theory the Hamiltonian was 
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taken as the positive relativistic energy of the free Weyl Hamiltonian. However, later it was shown 
that the adequate Hamiltonian should contain all the first order terms in Planck constant [8] which 
can be attained by employing the method introduced in [9]. Independently, the same Hamiltonian 
was conjectured in [TO] to restore the Lorentz invariance of the semiclassical chiral theory. To obtain 
it one first has to derive the Hamiltonian of the massive spin-1/2 fermion and then take the massless 
limit. 

Massive fermions also appear in condensed matter systems which were studied in terms of wave 
packets in mm- The semiclassical kinetic theory of Dirac particles was also discussed in m, where 
the Berry gauge fields have been given in a different basis and some classical degrees of freedom have 
been assigned to spin. We would like to establish the semiclassical kinetic theory of Dirac particles 
within the formalism given in (Tj. There are some advantages of employing this method. First of all 
because of not attributing any classical variables to spin but taking them into account by considering 
quantities matrix valued in “spin space”, the calculations can be done explicitly. The differential forms 
method provides us the solutions of the equations of motion for the phase space velocities in terms of 
phase space variables straightforwardly. Thus the particle currents can be readily derived. Moreover, 
we will show that within this formalism one can study the relativistic correction known as Thomas 
precession M- 

Thomas precession stems from the fact that a Lorentz boost can be written as two successive 
Lorentz boosts accompanied by a rotation which is called Thomas rotation. This purely kinematic 
phenomenon is essential to obtain the classical evolution of electrons spin correctly, without referring 
to the Dirac equation. Thomas precession should also contribute to the equations of motion of phase 
space variables. In fact, due to Thomas precession the covariant formalisms of Dirac particles yield 
equations of motion where anomalous velocity terms do not emerge m I16j . However, as we will 
see the equations of motion derived within the wave packet formalism possess anomalous velocity 
terms arising from the Berry curvature. This would have been expected because of the fact that our 
nonrelativistic formalism is not aware of Thomas rotation. Correction due to Thomas rotation should 
be installed in the formalism. We will show that our formalism suites well to take this correction into 
account: It contributes to the one-form obtained by the semiclassical wave packet on an equal footing 
with the Berry gauge field. In fact, it yields the cancellation of the anomalous velocity terms ignoring 
the higher order terms in momentum. The connection of Berry gauge fields and Thomas precession 
was first observed by Mathur m- 

In Section [2] the one-form corresponding to the first order Lagrangian is obtained by the wave 
packet composed of the positive energy plane wave solutions of the Dirac equation. Then the related 
symplectic form is constructed and the solutions of the equations of motion for the velocities of the 
Dirac particle are established in Section [3l Spin degrees of freedom are taken into account by letting 
the velocities be matrix valued. So that, one should consider a matrix valued distribution function. In 
contrary to to spin, helicity operator is a conserved quantity for the free Dirac particle. Moreover, we 
would like to split the particles as right-handed and left-handed appropriate to consider the massless 
limit and the chiral currents when there is an imbalance of chiral particles. Therefore we introduce 
a change of basis to the helicity basis as clarified in Section |4| Employing distribution matrix in the 
adequate basis we then can write the particle number density and the related current by the velocities 
written in terms of the phase space variables in Section [S] In Section [5] we obtained the massless limit 
by constructing the helicity eigenstates explicitly. In Section [6| a brief review of Thomas rotation is 
presented. Then, we presented how it appears in the one-form obtained by the semiclassical wave 
packet. We will see that up to higher order terms in momentum it contributes as the the Berry gauge 
held but with an opposite sign. Our semiclassical formalism should be supported by an equation 
governing spin dynamics. This will be attained employing Gosselin-Berard-Mohrbach (GBM) method 
|9] which is also needed to derive the semiclassical Hamiltonian. In the last section the results obtained 
and some possible applications are discussed. Moreover, we clarihed in which frame we obtained the 
semiclassical theory. 
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2 Semiclassical wave packet 


Dirac particle interacting with the external electromagnetic fields 5, B, whose vector and scalar 
potentials are a{x) and ao(x), is described by the Dirac Hamiltonian H = Hq + eao{x), where 

Ho{p — ea{x)) = /3m + cx ■ {p — ea{x)). (1) 


We set the speed of light c = 1 and let e < 0 for electron. The representation of the ctj; i = 1,2,3, 
and 13 matrices are chosen as 


ai = 






where Uj are the Pauli spin matrices. 

A semiclassical formulation given in terms of the wave packet made of the positive energy solutions 
of the Dirac equation provides a well defined one particle interpretation. Let the position of the wave 
packet center in the coordinate space be Xc, and the corresponding momentum be Pc- The semiclassical 
wave packet is defined in terms of the positive energy solutions, u°‘{x,p)] a = 1,2, as 


^x{Pc,Xc) = '^CaU°‘{Pc,Xc)e 
a 

For simplicity we deal with constant coefficients. We would like to attain the one-form rj which is 
dehned through dS as 

dS= j [dx\5{x, - X)^i i-ihd - HoDClt) = 

a0 

Hqd is the block diagonal Hamiltonian which should be derived from ([T]). Now, calculating dS one 
attains the rj one-form as follows. 

= —6‘^^Xc ■ dpc — • dxc — • dpc — HQ^dt. ( 2 ) 


Here we introduced the matrix valued Berry gauge fields 

= ihu^^°‘\pc,Xc)-^u^^\pc,Xc), 

OXc 

= ihu^‘'‘^'^{Pc,Xc)^U^^\pc,Xc). 

Opc 

Although we deal with the (3 -|- l)-dimensional space-time, the derivation of rj does not depend on 
dimension. 


3 Semiclassical Hamiltonian dynamics 


Instead of solving the Dirac equation in the presence of the electromagnetic vector potential a{x), we 
substitute p ^ p + ea{x), in ([T|) and consider the free particle solutions with E = y'jp + rn?. Then 
the positive energy solutions will not possess x dependence. Therefore by renaming {xc,Pc) {x,p) 
and setting = 0, we obtain the one-form 

rj = Pidxi + eaidxi — Aidp^ — Hdt. (3) 


The repeated indices are summed over. We suppress the matrix indices a, (3, and do not write explicitly 
the unit matrix I, when it is not necessary. We deal with the Hamiltonian H = H^{p) + eao{x), where 


Hd{p) = E 


he 


m- 


<T B 

2^ 




{B-p){(t-p) \ 

2E‘^{E + m) ) ■ 


( 4 ) 
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This semiclassical Hamiltonian includes all contributions which are at the first order in h. It is attained 
by making use of the GBM method [9]. 

To obtain Hamiltonian dynamics we have to introduce a symplectic two-form. In general the 
constituents of the one-form (matrix) r] can be non-Abelian, so that we adopt the definition of the 
symplectic two-form Cj to be 


Cj = dr] — —7] A 77. 
h 


We would like to emphasize that it is a matrix in spin indices. Employing the one-form ([3]), it yields 


Co = dpj^ A dxi + eSidxi f\dt + fidp^ A dt — G + eF. 


Here, £i = — , is the electric held and 


We always keep the terms up to the hrst order in h. The two-forms F = ^Fijdxi A dxj and G = 
^Gijdpj^ A dpj are written in terms of the electromagnetic held strength 


and the Berry curvature 




^ijk^k j 


G 




dpi 


dAj 

dpj 


i[Ai,Aj\ 


(^ijkGk- 


( 6 ) 


In order to derive the Berry gauge helds let us present the positive energy solutions for Hq{p) = 
a ■ p + f3m. They can be obtained as 


u“(p) = C/(p)uo“, 


where = (1 0 0 0), = (010 0), are the rest frame solutions and U{p) is the Foldy- 

Wouthuysen transformation: 

U(P)= + 


y^2E{E + m) 


Now, the Berry gauge helds can be written as 


dpi 

/+ projects onto the positive energy subspace. Hence, we acquire 

(T X p 


A = -h 


By making use of it in ([UD one gets 


G = - 


hm 


2E{E + m) 

p{a ■ p) 


cr + 


m{m + E) J 


Observe that the covariant derivative of G vanishes: 


riG- i 

d,g, = ^-^-[A,g,] = 0. 


( 7 ) 


( 8 ) 


(9) 
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To derive the equations of motion let us introduce the vector field 


d ^ d d 

Ot OXi opi 


( 10 ) 


{xi,pi) are the matrix-valued time evolutions of the phase space variables (xi,pi). By demanding that 
the interior product of the vector field with the symplectic form satisfy 


ifjUJ = 0, 

the equations of motion are found to be 

Xi — /i T Pj Gjii 

Pi = eSi - e Xj Fji. 

Time evolution of the volume form 

= -^oj A oj A 00 A dt, 

will yield the Liouville equation. The volume form can be written in terms of the canonical volume 
element of the phase space dV as 

Q = Coi/2 dV A dt. 

Here Wi /2 is the Pfafhan of the following (6 x 6) matrix, 




Obviously, 00 - 1/2 is still a matrix in the spin indices. We will present its explicit form. 

The Liouville equation can be provided by calculating the Lie derivative of the volume form 
associated with the vector field (HOI) . It can be calculated in two different ways. First, the Lie 
derivative of the volume form can be written in terms of the Pfaffian as 


= {iyd + diij){oo-i/2dV Adt) 


ddo^l2 ^ d[xiGo^i2) ^ d{Cb^,2Pi) ' ^ 


dt 


dxi 


dpi 


( 11 ) 


The second way is 


1 o 

C/jQ. = [iyd + div){—oo'^ A dt) 

- 


( 12 ) 


Now, calculate (I12p explicitly and compare with (lllh . This provides the Pfaffian as well as the velocities 
of the phase space variables in terms of the phase space variables themselves: 


hm 


dOi/2 



„ hem ( „ 

XOOi/2 



„ „ hAm 

W 1 / 2 P 

= e£ — ef X B - — 

■’ 2E^ 


m{m -|- E) 

{£ X p){(T-p) 

m{m -|- E) 


-eB{f-G), 


mym -|- E) 


On the other hand the term which survives in the Lie derivative is 

e^m{a ■ p){S ■ B) 


= h- 


E^{E + m) 


-dV A dt. 


(13) 

(14) 

(15) 

(16) 


This can be observed by explicit calculations. 
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4 Helicity basis and the continuity equation 


consider a matrix valued distribution function. For a Dirac particle spin is given by S = 


We study spin-1/2 massive particles where the velocities are matrices in spin indices. So that we should 

O' 

0 cr^ 

Hence, for the semiclassical wave packet composed of the positive energy solutions spin matrices 
become 51 = <t“^. They do not commute with the free Dirac Hamiltonian, so that they are not 

conserved in time. However the helicity operator 




u^ = 


cr ■ p 


P P 

commutes with the free Dirac Hamiltonian. It is appropriate to split up the particles as right-handed 
and left-handed which can be performed in the basis where the helicity is diagonal. In order to 
establish the diagonal basis we use the spherical coordinates where A is diagonalized by the unitary 
matrix 


R = 


cos(|) 

sin(f)e*v’ 


-sin(|)e-*^ 

cos(|) 


Thus the helicity basis is defined by 


(j) = Ru. 

Now, we can dehne the distribution matrix in the helicity basis as 


/0 = 


/r 0 
0 fL 


Inverse transformation will lead to the distribution matrix in the initial basis: 

-^2 e , f ■„2 e 


fu = RUR^ = 


/ijcos^ I/l sin^ I IsinOe fi) 

IsinOe'-^ifR - Jl) /rcos^ |-h /l sin^ | 


Since, for massive spin-1/2 particles the number of the right-handed and left-handed particles are 
equal we set /r = /l = /• Therefore, in the initial basis the distribution matrix becomes 

fu = /I. 

Let the distribution function, /, satisfy the collisionless Boltzmann equation: 

df 


dt 


df df ■ df ■ 

dt dxi dpi 


(17) 


In order to write the continuity equation, one should identify the particle number density n{x,p, t), 
and the particle current density j{x,p, t). However, we should first give an appropriate definition of the 
classical limit. This is done by taking the trace over spin indices and defining the classical velocities 
as follows 

\lW = Tr [wi/ 2 ], a/it a; = Tr [iwi/a], VWp = Tr [wi/ap]. 

We can write the probability density function as p{x,p,t) = VWf. Hence, the particle number density 
and the particle current density are given by 


n{x,t) = 


/ 


df’p 

(27r)’ 


Tr[wi/2]/, j{xR) = 


d^p 

(27r; 


gTr [iwi/ 2 ]/. 


Employing m and dH]) and setting f j 2 ^ —Op 
contribution from the the momentum space boundary, we find 


dp d(uji/ 2 Pt) _ because we suppose that there is no 


dn . 


/ 


= 0, 


(27r 


E^{E + m) 
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Thus we established the continuity equation for Dirac particles. 
One can also obtain the particle current 


j = 


d^p 

(27r)3 


Tr [iwi/z]/ 


d^p p 


5 Massless Fermions 


The massless Dirac equation can be written as two copies of the Weyl equation one for right-handed 
and one for left-handed fermions. Therefore, the helicity basis is suitable to discuss the massless case. 
First of all in the limit of vanishing mass, the Berry curvature in the initial basis ([SD yields 

Go = -hb^^. (18) 

P 

We introduced b = pj2p^ which is the monopole held situated at the origin of momentum space: 
V • 6 = 27rd^ {p). As it has already been announced the helicity basis is appropriate to discuss the 
massless case. In fact, when we change the basis from u to </>, the Berry curvature (I18p becomes 

Go^ = R^ GoR = -hb(^^ 

Now, the Berry curvature is singular and instead of ([9|) it satisfies 

= - 277 ha J^{p). 
dpi 

One can deal with the right-handed and left-handed fermions independently. Let us consider the 
right-handed massless fermions in terms of the projection 

=(i o) • 

Projection of the Pfaffian, Pr ojPr yields the scalar value 

y/Zj = 1 — ehb ■ B. 


Projecting (I14p and (11511 one acquires the solutions for the Weyl fermion as 

^/iox = —f^ — he£xb + ehB{f^-b), 

y/Zjp = e£ — ef^B + e^hb{E ■ B). 


denotes the massless limit of /, ([5]), projected by ([T9]). 

Let / be the distribution function of the right-handed fermions which satisfies the collisionless 
Boltzmann equation (|17p . The particle current density j can be written as 

^ ^ -he£xb + heB{r ' b)) f. 

The last term where the current is parallel to the magnetic field is known as the chiral magnetic effect 
term. The continuity equation becomes anomalous: 


, V . 1 

dt ^ ^ J 



0,t)£ • B. 


The Berry monopole situated at |p| = 0 is responsible for the non-conservation of the chiral particle 
current. 
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6 Thomas Precession 


We would like to clarify how to take into account Thomas precession [T3] within the wave packet 
formalism as presented here. To this end let us first briefly recall Thomas precession following [18j . 
The source of this phenomenon lies in the fact that if one would like to write a Lorentz boost as 
two successive Lorentz boosts, she should also rotate the coordinates with an angle depending on the 
velocities. This rotation yields an angular velocity known as Thomas precession. 

Suppose that a particle is moving with the velocity v with respect to laboratory frame at time t. 
Hence the particle’s co-moving frame denoted by the inertial spacetime coordinates x', is connected 
to the spacetime coordinates of the laboratory frame x, by the Lorentz boost Xboost{v) at time t : 

X Xboosti'^}^' 

Let the particle accelerates, so that it moves with the velocity v -h dv with respect to the laboratory 
frame at time t-hdt. Then at time t + dt the co-moving coordinate frame coordinates denoted x", will 
be connected to the laboratory frame x, by the Lorentz transformation 

X Xboosti,'^ “b dv'jx. (^l) 

Now, let us write the connection between the two co-moving frame coordinates x' and x" as 

x" = Xtx' . 


Inspecting ()20p and ()21ll the transformation Xt can be written as 

Xt Xjyoosti,'^ “b dv'j Xboosti ^)* 

The Lorentz boost Xboosti'i’ + dv) can be separated into two successive Lorentz boosts accompanied 
by the rotation R{dO) 

Xboosti,'^ “b dv'j — Ri^dO'j Xboost(,dv) Xboosti,'^)' 

The infinitesimal angle of rotation is 

dO =- V X dv. (24) 

7-1-1 

7 is the Lorentz factor. Therefore by plugging (I23p into (I22p one obtains 

Xt = R{dO)Xboost{dv). (25) 

When one deals with nonrelativistic equations of motion it is desired that in contrary to x”, 
successive co-moving frames of the particle would be connected only by boosts without any rotation. 
Therefore, x'" coordinates of the frame moving with the velocity v + dv, at time t + dt, will be obtained 
from the system moving with the velocity v, at time t, only with the boost Xboost{dv) without any 
rotation: 

X Xboostid'^')X . 

Then by making use of ()20p and ()23p the coordinates of the co-moving frame, x"', are written in terms 
of the laboratory frame coordinates as 

x'" = R{-dG) Xboost{v + dv) X. (26) 


How this phenomena manifests itself in the semiclassical kinetic theory of Dirac particles? In the 
wave packet formalism one deals with the group velocity 


_ dE p 
dp E 


(27) 
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Since 'y = Efm, hy plugging (1271) into one gets 

m{E + m) 

Let the laboratory frame and co-moving reference frames coincide at the time t = 0 when the particle 
is at rest. Hence the solution of the Dirac equation in laboratory frame at t = 0, is rt(0). Now, in the 
nonrelativistic kinetic theory formulation one deals with 


Mp) nr = u"'{p + dp) - u'{p) 

— ./?( do) ^boostip E dv)u(0) ^boost('^)R(0) 

= R(—dO)u(p + dp) — u(p). 

Here we ignore the terms at order of so that we write dv = dp/m. Therefore, when one considers 
Thomas precession, the Berry gauge field related term in the r/ one-form ([2]), of the nonrelativistic 
formalism will be 

ihv/{p)du{p)^j^ = ihv/{p) [R{—dO)u{p + dp) — u{p)] (28) 

The infinitesimal rotation of spinors is given by 

R{dO) = 1 - ^aijduj^^ = 1 + V{de) 

We implicitly work with the metric g = diag(l, —1,—1,—1), so that we have 

du^^ = R^^d9^ = -e^^^d0^. (29) 

Therefore, the rotation is expressed by 


vm = / 


" ). 

a ■ dOj 


Then, keeping only the first order terms in dp, (I28|) yields 

ihu\p)du{p)^j^ = ihv/■ dp + ihv){p)V{—d6)u{p). 

The first term is the Berry gauge field calculated in ([7]). The Thomas correction term can be shown 
to be 

ihu\p)V{—dO)u{p) = ■ dp. (30) 

4m ^ 

Ignoring the terms at the order of p‘^ the Berry gauge fields ([7|) yield the same contribution as in (|30p , 
up to a minus sign. Therefore, when the Thomas correction is considered the velocities can be read 
from (IT^ . (fTTp and (IT^ by setting the Berry gauge fields to zero. We conclude that the anomalous 
velocity terms disappeared when Thomas rotation is taken into account. This is in accord with the 
results obtained in relativistic formulations of Dirac particles MM- 


7 Time Evolution of Spin 

In the semiclassical approach adopted in this work, there is no classical degrees of freedom correspond¬ 
ing to spin of the particle. However, for the sake of completeness we would like to supply the time 
evolution equation of spin matrices. In Section U] we found that for the wave packet composed of the 
positive energy solutions spin is given by the Pauli spin matrices, cr. 
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To furnish the semiclassical equation governing time evolution of spin we may employ the GBM 
method [9]. In the GBM method one deals with the semiclassical phase space coordinates which are 
noncommuting: 

axp 

r = X — n——— -r. 

2E{E + m) 

Now, we can define the time evolution of spin by 



da 

dt 


ih 


[a,HD{p) + eao(r)] 




e[a,ri] 


dao(x) 

dxi 


(32) 


We deal with time independent electromagnetic potentials, so that {dao/dx) = S. Moreover, in the 
Hamiltonian Q we ignore the higher order terms in : 

Hd{p) = E-^-^^. (33) 

By plugging (f3T]l and (l33t) . into (1321) and setting j = E/m, one establishes the time evolution of spin: 


da e 
—— = —a X 
dt m 


-B + ^ 

7 7 + 1 


E X V 


(34) 


Observe that 


is the Bargmann-Michel-Telegdi equation m as composed in [18 


8 Discussion 

Semiclassical kinetic theory of massive spin-1/2 particles is studied within the method where the 
symplectic two-form is a matrix in the so called spin indices which label the positive energy solutions 
of the Dirac equation. These solutions are the building blocks of the wave packet which leads to the 
semiclassical approximation. Differential forms are employed to define the semiclassical Hamiltonian 
dynamics of Dirac particles. Solutions of the equations of motion for phase space velocities in terms 
of phase space variables are derived. To define the continuity equation of the particle number density 
and particle number current one has to define the distribution matrix adequately. This is possible in 
the basis where the helicity operator is diagonal. Therefore, a change of basis is performed. This is 
also needed to establish the massless limit. In the massless limit the expected anomalous continuity 
equation as well as the particle current possessing the chiral magnetic effect are obtained. 

Thomas precession correction needed in the nonrelativistic formulation of dynamics is studied 
within the wave packet formalism. We showed that up to higher order terms in momentum it sweeps 
out the contributions arising from the Berry gauge fields. This coincides with the results obtained in 
relativistic formulation of Dirac particles. This is the main result obtained in this work. The presented 
method of introducing Thomas precession correction is valid in general. It can be applied to other 
semiclassical approaches of Dirac like systems where the underlying Hamiltonian of the theory is given 
by Dirac like Hamiltonian as in some condensed matter systems. 

We built the semiclassical kinetic theory of Dirac particles starting with the wave packet. This 
formalism is valid when the potentials do not vary rapidly across the wave packet. Now, the particle 
is an extended object whose diameter is equal to the Gompton wavelength [12]. It is a self-rotating 
wave packet, so that its center of mass and center of charge or the “center of energy,” do not coincide 
[20] . However, the circular current which yields this interpretation is derived from the equations of 
motion for the center of mass P. Thus, in the absence of the electromagnetic vector field a = 0, 
kinetic part of the one-form (l3|) can be written in the symmetric form: 

7k = 2 (PcM ■ dxQyi — ®CM ■ dpcm) ~ A ■ dpcM- (35) 


10 








As we have already observed in Section [U spin of the wave packet is given by S' = |<t. Thus the Berry 
gauge field can be expressed as 


A = —— - S X V ^ - S X V, 

E + m 2m 

by keeping the terms linear in the velocity v. On the other hand for a spinning object the center of 
mass and the center of energy are related as 


®CM ~ ®CB H S XV — iCcB ~ 2A, 
m 

PcM ~ PCE- 


(36) 

(37) 


Therefore, in terms of the center of energy variables (xce,Pce), tbe kinetic one-form (j35p . can be 
rewritten in the canonical form: 

?7k ~ 2 (PcE ■ dxQE ~ Xqe ■ dpc^) ■ (38) 

The Thomas precession affects the center of mass and leaves intact the center of energy [20]. Phase 
space variables of the GBM method refers to the center of energy. In fact because of this one can 
directly take the vanishing mass limit [2T]. Moreover, as we showed in Section [T] the Hamiltonian 
density dH), which depends only on the momenta pcE ~ Pcm, engenders the correct time evolution of 
spin without adding a Thomas precession term. 

The semiclassical kinetic theory formulation of Dirac particles and obtaining the massless limit by 
constructing the suitable basis which we developed here, can be generalized to systems where some 
other interaction terms are present in the underlying Hamiltonian. 
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